Tunneling Anomaly in Superconductor above Paramagnetic Limit 
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We study the tunneling density of states (DoS) in the superconducting systems driven by Zeeman 
splitting Ez into the paramagnetic phase. We show that, even though the BCS gap disappears, 
superconducting fluctuations cause a strong DoS singularity in the vicinity of energies —E* for 
electrons polarized along the magnetic field and E* for the opposite polarization. The position of 

the singularity E* = | \Ez + \J E z — A 2 ^ (where A is BCS gap at Ez = 0) is universal. We found 

analytically the shape of the DoS for different dimensionality of the system. For ultrasmall grains 
the singularity has the form of the hard gap, while in higher dimensions it appears as a significant 
though finite dip. Our results are consistent with recent experiments in superconducting films. 
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It is known that magnetic field, H, suppresses super- 
conductivity because it violates the time reversal sym- 
metry |Q. Typically, this symmetry breaking is associ- 
ated with effect of the magnetic field on the orbital mo- 
tion of electrons. However, in some physical situations, 
the main mechanism of the destruction of the Cooper 
pairing is the Zeeman splitting of states with the same 
spatial wave-functions but opposite spin directions. One 
can consider, e.g., a thin superconducting film placed in 
magnetic field parallel to the plane of the film 0. Re- 
cently, another possibility was exploited experimentally 
H - ultrasmall superconducting grains. In both cases, 
the size of a Cooper pair is restricted geometrically and 
a flux through this pair reaches flux quantum at fields 
higher that the bulk critical filed H C2 . We assume that 
superconductivity is already destroyed by the Zeeman 
splitting Ez = 9l^bH when it happens jjj, (here gr, is 
the g-factor Lande, and (j,b is the Bohr magneton). 

Strictly speaking, the spin splitting destroys supercon- 
ductivity as soon as Q Ez > \/2A, with A being the su- 
perconducting gap. This transition from superconduct- 
ing to paramagnetic normal state is of the first order, and 
these phases coexist in the interval A < Ez < 2A. From 
now on, we assume the condition Ez > V%A. 

One might expect that the Cooper pairing is irrele- 
vant for the properties of the normal paramagnetic phase. 
In this Letter we show that, on the contrary, there are 
clear and observable effects of the pairing in paramag- 
netic state even far from the transition region. 

One of the most fundamental manifestation of the su- 
perconductivity is the gap in the tunneling density of 
states (DoS) around the zero energy ||. This gap ap- 
parently disappears when the system becomes paramag- 
netic. We will show that at the same time there appears 
a dip in the DoS. The shape of this dip depends on the 
dimensionality of the system. However, its position E* 
is remarkably universal: 



for OD (grain), ID (strip) and 2D (film) cases. 

This result should be compared with the the singu- 
larity in the DoS due to the usual superconducting fluc- 
tuations in the normal metal. According to Ref. Q, in 
the magnetic field another anomaly appears in addition 
to the zero bias anomaly when bias V corresponding to 
Zeeman splitting eV s = Ez- The singularity we consider 
here is positioned at substantially lower energy and, as 
we will see, is stronger than those considered in Ref. [Q. 

Recently, singularity of this type was observed in gran- 
ular Al film in a parallel magnetic field Q|. We believe 
that the deviation from the law eV s = Ez observed in 
Ref. is coherently explained by our Eq. (|l|). 

We begin with the simplest but instructive case of OD 
system (ultrasmall grain). The Hamiltonian TL of this 
system can be written in the basis of the exact single- 
electron states for non-interacting system as 
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Indices i,j and a =|, J, label the orbital and spin state of 
an electron: its total energy Ei^ a is the sum of orbital e, 
and spin parts, -E^-m) —^iTEz/2; and a ia ,Oi^ are the 
corresponding fcrmionic creation-annihilation operators. 
Finally, 6 = (e;+i — £i) is the mean level spacing and 
A is the dimensionless interaction constant. In Eq. (||), 
we omitted some diagonal terms (those with orbital in- 
dices equal pariwise), which with the help of Hartree- 
Fock approximation can be included into the definition 
of the eigenenergies ei. We also omitted off-diagonal 
terms (involving the fcrmionic operators with non-paired 
indices): the corresponding matrix elements are known 
Q to be smaller than the diagonal ones by a large fac- 
tor 1/g, where g 3> 1 is the dimensionless conductance 
of the grain. Hamiltonian (^) is nothing but the usual 
BCS Hamiltonian, and it was used in numerous publica- 
tions H on the properties of ultrasmall superconducting 
grains. 
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BCS instability which system (Q) has at temperature 
T = and H = 0, disappears as soon as Ez exceeds 
\/2A. In the absence of the superconducting gap, the 
structure of the ground state is similar to that without 
interaction: orbitals with energies > Ez/2 are spin- 
less (orbitals with a < —Ez/2 are double ocupied while 
those with e, > Ez/2 are empty), while orbitals in the 
energy strip |e»| < Ez/2 are spin polarized with spin up 
(t) (we measure all the energies from the Fermi level). 
The interaction term in the Hamiltonian (0) does not 
affect the spin polarized states, but mixes the double- 
ocupied and empty states. Those states are separated 
from each other by large gap Ez- Therefore, this mix- 
ing is perturbative, and it does not change ground state 
qualitatively. 

Contrarily, the spectrum of the excitations, e.g. the 
tunneling DoS changes drastically due to the interac- 
tion. Consider a spin-down electron with energy < 
E < Ez/2 entering the grain. The orbital energy of this 
electron e is close to eo = E — Ez/2, and this orbital is 
ocupied by a spin-up electron. Therefore, the tunneling 
event creates the spin-singlet state with the energy 2eo, 
and this state can mix with all the empty states ej. This 
mixing turns out to be resonant at some energy E = E*; 
for E close to E* it requires nonperturbative treatment. 

Before turning to the rigorous calculations, let us dis- 
cuss this effect qualitatively using the following simpli- 
fication. Instead of the whole many-body system, we 
consider two-electrons only, however, the single-electron 
orbitals for this pair are restricted by the orbitals with 
£i > Ez and by one orbital eo. The role of the rest of 
the electrons is to restrict the Hilbert space for a given 
electron pair. This simplification is similar to the Cooper 
procedure |10|. The interaction in Eq. (Q) involves spin 
singlet orbitals only. Thus, the wave function of the elec- 
tron pair ip can be labeled by one orbital index and it 
obeys the Schrodinger equation sipi = 2e^ — A<5 ^j- 
The eigenenergies e of this equation can be found from 
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For low-lying eignestates e < Ez one can substitute the 
summation in Eq. (j^) by integration. Given the high- 
energy cut-off Q, it yields 
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For small level spacing 5 <C Ez, A;,, the pair energy s is 

1/2 
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n b = E z -A b . (5) 



Prior the spin down electron tunnels in, the energy 
of the spin-up electron on the orbital eo was equal to 



Ei = eo — Ez/2, and thus the energy of one electron 
excitation E — e — E± is given by 



E = E* b ± 
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The origin of tunneling anomaly is now transparent from 
Eq. (^): due to the repulsion between the state formed 
immediately after an electron tunnels into the system 
and the bound states of the Cooper pair, there is no 
spin-down one electron excitations in the energy strip 
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hard gap in the DoS is formed. 



It is important to emphasize that (i) the width of this 
gap y/SSAf, significantly exceeds the single electron level 
spacing and (ii) this singularity persists even when the 
system is deep in the paramagnetic phase Ez S> \ 2A. 

It is also noteworthy that if a spin-up electron tun- 
nels into the grain, it never finds the pair for itself, and, 
therefore, no tunneling anomaly happens in this case. It 
means that the overall DoS does not vanish but rather 
shows the suppression by a factor of two. However, for 
the spin-polarized electrons tunneling into the grain, we 
predict the complete suppression of the tunneling DoS. 

The same arguments allow to justify the similar sin- 
gularity, when spin-up electron with energy —Ez/2 < 
E < tunnels out from the system, while the spin down 
electrons tunneling from the system are not affected. 

The qualitative consideration above grasps the correct 
physics, however it fails to describe the effect quantita- 
tively, it predicts correctly neither the position [compare 
E*, Eq. (§), with Eq. (Q)] nor the width of the gap. This 
is similar to the discrepancy between the binding energy 
A b in the original Cooper procedure and the correct BCS 
gap A: all the electrons below the Fermi energy were 
frozen. To remedy this drawback, we employ a paramet- 
rically exact procedure described now. 

We start with the propagator of the superconducting 
fluctuations. The diagrammatic equation, Fig. [l]a, yields 
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Here A = we ' is the BCS gap, and u+ 
Propagator (fjj) has the pole at ui = ±0: 



(7) 



bj + iOsgnw. 



n 



E z 



A 2 



(8) 



n has the meaning of the bound state energy of two quasi- 
particles. Using Eq. (||) instead of oversimplified Eq. (|5|) 
in the expression for E£ in Eq. (||) , we obtain Eq. (Q) for 
the position of the singularity E* . 

Consider now the shape of the tunneling DoS 
near the singularity. DoS for spin polarizations a =f 
, I can be expressed through the one particle Green 
function (GF) iGi >a (uj) at zero temperature: v a (u>) — 
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— l-sgnwlm^ Gi In its turn, the GF is given by 
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the one-particle self-energy 11 1. 

The leading part of Sj.j., diagram Fig. [j]b, is given by 
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The singular contribution to the integral in Eq. (^|) comes 
from the positive pole of A c , which corresponds exactly to 
the repulsion between the bound state of the Cooper pair 
with the state formed after tunneling which we discussed 
earlier. From Eq. ([)]) we obtain for the GF 
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Summing up over all the orbitals and neglecting the 
fine structure of the DoS on the scale of 8, we find 
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where vq is the bare density of states per one spin. Equa- 
tion (O) describes the hard gap of the width 
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positioned at w = _E* given by Eq. ([!]). 
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FIG. 1. Diagrams describing (a) the fluctuation propaga- 
tor A c (w) and (b) first order contribution to the self energy 
Si,|(w). Diagrams (c) and (d) show second order contribu- 
tions to E, which were neglected in comparison with reducible 
diagram (e). 

Higher order corrections to the self-energy (e.g. shown 
in Fig. |l]c,d) are negligible. Indeed, let us compare the 
contributions of those diagrams with the reducible dia- 
gram Fig. [l]e which is included in Eqs. (j|) and (|To|). At 
lo = E*,2ei = ft, all three GF in diagram (e) diverge 



after integration over intermediate frequencies, whereas 
only two GF diverge (non-resonant GF are indicated by 
• ) in Fig. 0c,d. It means that diagrams (c-d) are smaller 
than the main contribution by a factor 8/ A -C 1. 

According to Eq. (pi]), the width W of the hard gap 
is much larger than the mean level spacing 8. However, 
it apparently vanishes as 8 1 / 2 for infinite systems. To be 
more precise, our 0D consideration was valid only if the 
width of the gap Wo ~ VTK does not exceed the Thou- 
less energy Et — hD/L 2 ( D is the diffusion coefficient). 
Since 8 ~ l/L d , (d is the dimensionality of the sample) 
the condition Wo < Et breaks down for bulk samples, 
L — » oo, in all physical dimensions d < 4 and our zero- 
dimensional description (|^) becomes invalid. Neverthe- 
less, the strong singularity in the vicinity of E* persists 
in the tunneling DoS for d = 1, 2 systems [fl2|| . 

We describe the interaction by usual Hamiltonian 
"Hint = — t , _1 A J drajetjajaj. In the bulk system, the su- 
perconducting fluctuations can be inhomogenous. Thus, 
the propagator similar to Eq. (0) depends on the wavevec- 
tor Q. Diagram for such propagator averaged over disor- 
der is similar to Fig. [|a and standard calculation yields: 
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A c at small Q, (DQ 2 <C f2) has the singularity at u> close 
to the energy from Eq. (||). Therefore, one may expect 
the singularity in DoS at the same energy E* . 

Evaluation of the first order correction to the 

DoS, similar to Fig. [l]b, confirms this expectation. Inte- 
grating over intermediate frequency, averaging over dis- 
order, and assuming DQ 2 <C ft, we find the singular part 
of the correction [Of 
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where d = 1, 2 is the dimensionality of the sample, and 
the Cooperon C CT is given by C^n = [lu ± E* + iDQ 2 ] . 
For ID case (strip), integration in Eq. (|l4|) gives: 



vo 



A 1 ' 2 



2v ttD 1 /2 \2\uj±E* 



A 



3/2 



(15) 



In two dimensions the correction to DoS vanishes in 
the lowest order. Second order correction similar to that 
in Fig. |l|e averaged over disorder reads [|l3) 
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For two dimensional film integration yields 
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where g = AnhDvo 3> 1 is the dimensionless sheet con- 
ductance of the film in the normal state, ft is worth 
noticing, that the singlularity given by Eq. ([if]) is much 
more pronounced than that due to the superconducting 
fluctuations in the normal metal [Q, oc g~ x hi[hi(u;±Ez)]- 
Both corrections ([Hi]) and ( [l6| ) diverge with approach- 
ing the singularity point to — > ±E*. It means that the 
higher order terms have to be summed up. Here we 
present only the results of this nonperturbative treat- 
ment |0| . The resulting DoS can be written as 



TU) 



1,2, 



(17) 



where the characteristic widths of the singularity W\_2 in 
one and two dimensons are given by 
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In one dimension DoS acquires a universal shape: 
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where y(x) = \/l + ix 3 . In two dimensions shape of the 
DoS weakly depends on the conductance g 

1 - z{x) 1 



3-2{x) = Rc- 



(19b) 



l + z(x)' — ix + ln(4gz) 

At x » 1, Eqs. ( |l9| ) reproduce perturbative results (JIB] ) 
and (H). DoS near the energy E* is shown in Fig. B. 
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FIG. 2. Singularity in DoS for spin-down polarized elec- 
trons for (a) ID, and (b) 2D systems. Widths of the singu- 
larity Wi,2 are given by Eq. (181), and the shape is defined by 
Eqs. (O) and M). 



Finally, let us mention that the DoS singularity is su- 
pressed by a finite temperature and by the spin-orbit 
scattering: the singularity remains observable only if 
W > T,h/r so . Corresponding expressions will be pre- 
sented elsewhere p3| . 

In conclusion, we have shown that the tunneling DoS 
of superconductors above the paramagnetic limit has a 
singularity at energy determined by Zeeman splitting and 
superconducting gap, Eq. (|l|). This phenomenon is not 
take into account by usual BCS mean field and can not be 
obtained within Gorkov-Nambu formalism. The shape of 
the singularity, obtained nonperturbatively, depends on 
the dimensionality of the system, Eqs. (|l7| ) and (|l9|). We 
believe that our theory is adequate for existing experi- 
ments H, and suggests new effect in ultrasmall grains. 

Discussions with A.I. Larkin and B.Z. Spivak are ac- 
knowledged with gratitude. 
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